Beware the two-beam coupling phase shift: Interpreting pump-probe measurements of 

the nonlinear refractive index 

J. K. Wahlstrand, 1 J. H. Odhner, 2 ' 3 E. T. McCole, 2 * 3 Y.-H. 
Cheng, 1 J. P. Palastro, 1 R. J. Levis, 2 ' 3 and H. M. Milchberg 1 

1 Institute for Research in Electronics and Applied Physics, 
University of Maryland, College Park, MD 20742, USA 
2 Department of Chemistry, Temple University, Philadelphia, PA 19122, USA 
''Center for Advanced Photonics Research, Temple University, Philadelphia, PA 19122, USA 

Nonlinear optics experiments measuring phase shifts induced in a weak probe pulse by a strong 
pump pulse must account for coherent effects that only occur when the pump and probe pulses are 
temporally overlapped, ft is well known that a weak probe beam experiences a greater phase shift 
from a strong pump beam than the pump beam induces on itself. The physical mechanism behind 
the enhanced phase shift is diffraction of pump light into the probe direction by a nonlinear refractive 
index grating produced by interference between the two beams. For an instantaneous third-order 
response, the effect of the grating is to simply double the probe phase shift, but when delayed 
nonlinearities are considered, the effect is more complex. A comprehensive treatment is given for 
both degenerate and nondegenerate pump-probe experiments in noble and diatomic gases. Results 
of numerical calculations are compared to a recent transient birefringence measurement [Loriot et 
at, Opt. Express 17, f3429 (2009)] and a recent spectral inteferometry experiment [Wahlstrand 
et at, Phys. Rev. A 85, 043820 (2012)]. We also present results of two new experiments using 
spectrally-resolved transient birefringence using 800 nm pulses in Ar and air and degenerate chirped 
pulse spectral interferometry in Ar. Both experiments support the interpretation of the negative 
birefringence at high intensity as arising from a plasma grating. 



I. INTRODUCTION 

The intensity-dependent refractive index due to odd- 
order optical nonlinearities is a fundamental phenomenon 
in nonlinear optics. Accurate knowledge of it is impor- 
tant for most applications of intense laser pulses. In con- 
densed matter, the Z scan technique, a measurement of 
self-focusing [1], is widely employed, but in gases, it is dif- 
ficult to use because the interaction length cannot be eas- 
ily characterized. Experiments that measure the phase 
shift induced by a pump pulse on itself in an extended 
medium [2H5] also require careful consideration of prop- 
agation effects and plasma defocusing at high intensity 
[BJ. Accounting for these phenomena quantitatively re- 
quires extensive theoretical modeling. For this reason, 
pump-probe experiments, which measure the effect of 
cross phase modulation, have an important advantage. 
For a weak probe one can at least be certain that the 
response of the medium is linear in the probe field, and 
interpretation of the experimental results appears rela- 
tively straightforward. However, in such an experiment, 
a pump-probe interference grating always appears in the 
medium and pump light is diffracted by this grating into 
the probe beam direction. One must properly account 
for the contribution of this to the probe phase shift. 

Depending on the relative phase between the diffracted 
pump beam and the probe beam, the probe can experi- 
ence a phase shift and/or an amplitude change [7]. This 
is the physical mechanism responsible for the enhanced 
phase shift in cross phase modulation compared to self- 
phase modulation. The amplitude effect is often referred 
to as two-beam coupling [8]. For an instantaneous optical 
Kerr nonlinearity, one finds that the probe picks up twice 



the nonlinear phase shift of the pump, but there is no net 
energy transfer. In general the result depends on the de- 
tails of the nonlinearity, such as the response time, and 
even the details of the pulse shape, such as whether or 
not it is chirped. A major goal of this paper is to provide 
a comprehensive discussion of the phase shift imparted 
on the probe pulse and the underlying physics. 

We approach this topic in the context of a debate over 
the higher-order Kerr effect, ignited by a recent degener- 
ate pump-probe experiment in Ar, N 2 , and 2 reported 
by Loriot et al. [9j. A negative birefringence was ob- 
served at high pump intensity when the pump and probe 
pulses overlapped in time. It was concluded that the 
expansion of the nonlinear refractive index in powers of 
the optical intensity I, An = ri2l + n^I 2 + jiqI 3 + is 
dominated by higher-order negative terms at high inten- 
sity. Initially, the enhancement of the phase shift in cross 
phase modulation was not considered, and the Kerr co- 
efficients for the effect of a pulse on itself were corrected 
in an erratum 10J. However, this correction accounts 
only for the nearly instantaneous bound electronic com- 
ponent of the nonlinearity. An alternative explanation 
for the negative signal observed by Loriot et al. [9] can 
be found by considering the plasma grating diffraction 
effects when the pump intensity is high enough to ionize 
the gas [llj . Recently an experiment in air using 400 nm 
pulses [T2] found that the intensity dependence of the 
negative signal was consistent with a calculation of the 
plasma grating-induced birefringence based on multipho- 
ton ionization |llj . 

Here, we generalize the theory to allow use of any ion- 
ization model, account for grating effects from the molec- 
ular alignment component of the nonlinearity, and extend 
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the theory to handle pulses of arbitrary frequency and 
chirp. We also present new experimental results from 
degenerate pump-probe experiments demonstrating the 
probe phase shift originating from plasma grating diffrac- 
tion effects. We note that many recent experiments have 
used a nondegenerate probe [T2TfT8] . for which there is 
no plasma or rotational grating contribution to the non- 
linear phase shift [TTlll2j . as detailed below. We first de- 
scribe a theoretical framework for the calculation of the 
outgoing probe field in a general pump-probe experiment. 
We then numerically calculate the signal and compare 
the calculations to the Loriot et al. experiment [5], and 
to nondegenerate spectral interferometry measurements 
[131 HH [T6HT8] . Finally, we present new experimental re- 
sults using spectrally-resolved transient birefringence and 
spectral interferometry with a degenerate probe pulse. 



which is in general not linearly polarized, is decomposed 
according to A p = xA x + yA y . To allow for chirped 
pulses to be handled, A e and A™ c are complex. We ne- 
glect diffraction and nonlinear propagation effects of the 
pump pulse. The primary pump depletion mechanisms 
are ionization and pumping of rotational states, which 
are negligible for the intensities and interaction lengths 
used in the experiments to which we compare our calcu- 
lations. 

The interaction between pump and probe beams oc- 
curs near the waists of nearly Gaussian beam profiles, 
and in practical experimental implementations, the beam 
waist is much larger than the wavelength of the light. 
The transverse derivatives d 2 /dx 2 and d 2 /dy 2 in V 2 are 
thus negligible in the interaction region, and the left side 
of Eq. (JlJ , keeping only terms propagating in the probe 
direction, becomes (omitting the complex conjugate) 



II. THEORY 

Initially reported, as far as we are aware, by Chiao, 
Kelley, and Garmire |19j . the grating-induced probe 
phase shift ( "weak wave retardation" ) is related to "self- 
diffraction" effects that have been often encountered 
in experiments in condensed matter (for example, |20l - 
123] ). Two-beam coupling, which typically refers to en- 
ergy exchange between the pump and probe beams, has 
been studied in detail [24H28] , However, the probe 
phase shift does not seem to have been discussed except 
for instantaneous (722) and photorefractive nonlinearities 

[221 ESI- 

Our general approach is similar to previous calcu- 
lations of two-beam energy coupling [7 . We assume 
two incident optical fields, a pump ("excitation") pulse 
E e (r,i) = (l/2)A e (r,t)e lk '= r - i ^ t + c.c, and a probe 
pulse Ejf c (r,f) = (l/2)Aj? c (r, i) e i V r -i'M + c .c. In the 
experiments in gases that are the primary focus of this 
paper, the lowest optical resonances are at roughly 10 
eV, so the linear response of the medium is nearly in- 
stantaneous, and the effect of dispersion over a typical 
interaction length of a few mm is negligible. We include 
the linear response of the neutral, isotropic gas in a con- 
stant real refractive index no and group the rest of the 
response in a nonlinear polarization Pnl- We solve the 
wave equation 
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with the incident fields as an initial condition. 

In experiments, the probe beam is isolated from the 
pump beam spatially [SJII3 and/or by frequency filter- 
ing [T3TfT7] . In the calculation we are therefore inter- 
ested only in terms that propagate in the k p direction 
and oscillate at frequencies near uj p . This outgoing probe 
field is E p (r,<) = (l/2)A p (r, t)e' w " ■ r ' l ^ t + c.c. We de- 
fine kp — kpZ^ kg - — - cos (9 + ysin#), and linearly 
polarized incident pump and probe envelopes A e = xA e 
and A p nc = xA" lc + yA 1 y nc . The outgoing probe envelope, 
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The pulses simulated are many cycles long (40 fs 
and above), so the slowly- varying envelope approxima- 
tion (SVEA) may be employed, using \d 2 A p /dt 2 \ <C 
\uj p dAp/dt\ and \d 2 A p /dz 2 \ <C \k p dA p /dz\. Using k p = 
nouip/c, we have 
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where Pnl.p includes only terms in Pnl containing 
gikp-r-iwpt ^ ag Qn jy they significantly affect the probe 
field). 

We next derive expressions for the polarization Pnl,p 
due to the nonlinear interaction of the pump and probe 
pulses with the medium. We assume Pnl = Pol + Ppi + 
P ro t, where P e i is the electronic nonlinearity, P p i the 
plasma (free electron) nonlinearity, and P rot the nonlin- 
earity due to the rotational response (molecular align- 
ment) . 



A. Electronic nonlinearity 

Cross phase modulation of a weak beam due to a third- 
order nonlinearity is well known and is treated in nonlin- 
ear optics textbooks (e.g., [5])- We provide details here 
in order to make a clear connection to the more compli- 
cated nonlinearities arising from ionization and molecular 
alignment considered later. Assuming only a third-order 
contribution to the polarization and a uniform medium, 
we have 

(Pel)i(M) = JJJ X\jll(t-tl,t-t 2 ,t-t 3 ) 

x E J (r,t 1 )E k (r,t 2 )E l (r,t 3 )dt 1 dt 2 dt 3 , 
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where x is the third-order nonlinear response function. 
As all electronic resonances are far higher in energy than 
huj e and Ti<jj p , to a good approximation x (t — t±,t — 
t2,t-t 3 ) w x [3) S(t - t^Sit - t 2 )S(t - t 3 ). The nonlinear 
electronic component of the polarization then simplifies 
to 



(PciMM) 



X %E ] (r,t)E k (r,t)E l (r,t). 



(3) 



We use the total field E P 



En 



-, p in Eq. (131). We em- 
phasize that only the terms propagating in the k p direc- 
tion, which are proportional to e lkp ' r ""' 1,i , contribute to 
the signal in the experiment. The terms proportional to 
gike-r-jw,.* j eac j ^ Q cnan g es m t ne phase and amplitude of 

the pump field, which we neglect. The terms proportional 

to e^ke-kp)- 1 "-^^-"?)* ; an d e i(2k p -k,.)T-i(wp-w,,)t j ea( j 

to optical fields that do not reach the detector in properly 
designed experiments. In degenerate pump-probe exper- 
iments the probe beam is isolated spatially, and in non- 
degenerate experiments, where typically the pump and 
probe beams are collinear, the probe beam is isolated 
spectrally. In the SVEA, it is shown in Appendix A that, 
defining the Kerr coefficient n 2 — l2ir 2 Xx 3 xxx/{rioc), we 
have for the nonlinear polarization term (x component), 



47T d 2 (P cl ) x _ 2n n 2 u} p 
c 2 dt 2 c 2 



—tUpI e A x — 2i—^-A x 
dA x i a , B . r r , 



(4) 



where I e = n c/(87r)|A e | 2 is the pump intensity. Off- 
diagonal elements of % are responsible for the y com- 
ponent of the polarization, and in an isotropic medium, 
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Details are given in Appendix A. 
Putting Eq. Q in Eq. §2§, we have 

dA x n dA x 2n 2 ( . dl e \ 2n 2 dA x 

ih ~~df = ~ r p/e -^) A *-— Ie ^r> 

where we have used uj p /k p — c/hq. Defining A x — 
|^4a;|e I * x , one can easily show that the nonlinear change 
in phase as the probe propagates through the interaction 
region is 
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Thus, the effective Kerr nonlinear refractive index for the 
probe is found to be Au™ ss = 2n 2 I e , twice as large as the 
change in nonlinear index calculated for the pump beam 
on itself, A% = n 2 I e . The phase shift of the probe is en- 
hanced by a factor of two regardless of whether the probe 



is degenerate or nondegenerate with the pump. Similar 
calculations using higher-order response functions predict 
enhanced probe phase shifts for the higher-order Kerr ef- 
fect due to x (5) , X (7) , etc. [TO]. 

The factor of two enhancement of the Kerr phase shift 
of a weak probe compared with the phase shift of the 
pump due to self phase modulation can be explained 
physically by diffraction from a nonlinear index grating 
0[TT]. To understand this, it is helpful to consider the 
problem in a different, less rigorous way. The total in- 
tensity is, in the limit of a weak probe, 



I {r ,t) = 7 ^[\A e \ 2 + (A:A x 



zAk r— iAujt 



+ c.c.)], (7) 



where Ak = k. p — k e and Auj = u> p — oj e . Interference 
of the pump and probe beams causes a sinusoidal modu- 
lation of the total intensity (the terms in parenthesis) 
where the beams cross. When inserted into Au^- = 
n 2 I, we find An K = An s K + (An 9 K e Ak - r ' iAut + c.c), 
where An s K = n 2 I e is the "smooth" nonlinear index and 
An 9 K = n 2 noc/ (8tt)AIA x is the "grating" nonlinear in- 
dex. The polarization due to this nonlinear refractive 
index is P c i = [An^no/(27r)]E. The terms that contain 
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Taking the second derivative with respect to time and 
applying the SVEA leads to Eq. Q. The smooth com- 
ponent An s K produces the same Kerr phase shift in the 
probe field as the pump induces on itself, and the grating 
component An 9 K combines with the pump field to pro- 
duce an identical phase shift in the probe field, doubling 
the overall response. 

Pump light is diffracted by the nonlinear index grating 
into the probe path with just the right phase and am- 
plitude to produce a doubling in the probe phase shift. 
Probe light is also diffracted into the pump beam direc- 
tion, and no net energy transfer occurs. The extra probe 
phase shift is independent of the magnitude of the probe 
field, since the index modulation scales with the probe 
field. When uj e ^ u) v , the Kerr grating moves in time as 
well as space. The diffracted pump light at frequency uj e 
is shifted to exactly match the probe light at frequency 
ui p , and the factor of two enhancement in the phase shift 
still occurs. Note that the crossing angle 9 does not ap- 
pear in Eqs. ( 4|5 1; for uj e ^ uj p the same factor of two 
enhancement in the probe phase shift occurs even for 
9 = 0, where there is no spatial grating. 
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B. Plasma nonlinearity 

Ionization of atoms and molecules by the intense op- 
tical field frees electrons, causing the refractive index to 
change with time. We calculated cross phase modulation 
effects from this plasma response previously assuming a 
multiphoton ionization model . Here we generalize to 
the case of an ionization rate w(I) that depends only on 
the cycle averaged intensity /, allowing coverage of the 
multiphoton and tunneling limits [30] [3T] , as well as the 
intermediate intensity regime where laser filamentation 
typically occurs in gases. 

The presence of free electrons leads to a polarization 
term 



d 2 P p 

dt 2 



-JV e (r,i)E(r,t), 



where N e is the electron density. The ionization rate 
is a function of time and space, and yields the electron 
density through dN e /dt = N a w (/ (r, t)), where N a is the 
number of atoms or molecules per cm 3 . Using the total 
intensity given by Eq. Q we have 



w(I(v,t))=w(I e (v,t)) 
dw 



n c 
8tt 
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+ C.C. 



Ie(r,t) 



As in the previous section, the total free electron den- 
sity can be split into smooth and grating terms N e = 
N* + (7V|e iAk ' r + c.c), where the smooth term is 



N!(r,t) = N a l w(I e {r, t'))dt' , 
and the grating term is 



Ni(r,t) 
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dw 
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A* e (r,t')A x (r,t')e- iAut 'dt'. (9) 



7 e (r,t') 



Note that we have assumed the weak ionization limit 
(i.e. neglecting depletion of the neutral atom popula- 
tion) , which is a good approximation for the experiments 
we compare to. The integral in Eq. (JoJ) is suppressed 
when uj p ^ u) e . 

As described previously, we are interested only in the 



terms containing e l 



as only they contribute to the 



signal in the probe beam direction. These are 
d 2 P n i 



dt 2 
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(10) 



The smooth term causes a phase shift in the probe pulse, 
consistent with a change in index An = —N e /2N CT , 
where N CI is the critical density. The pump beam com- 
bines with the grating term Njj to create a polarization 



source propagating in the probe direction. In the calcu- 
lations shown later we shall find that the main result of 
the grating term is a negative phase shift for the probe 
polarization component parallel to the pump polarization 
[llj . The grating phase shift only appears when uj p ps u) e . 



C. Rotational nonlinearity 

In molecular gases, an intense laser field applies a 
torque to the molecules, partially aligning them with the 
optical field [32-35 . The effective third-order nonlin- 
earity due to this alignment can be an important, even 
dominant part of the total nonlinearity for intensities be- 
low the ionization threshold [HJ [XT] [35]. The nonlin- 
ear rotational response comes about because in diatomic 
molecules the polarizability depends on the angle be- 
tween the laser field direction and the molecular axis. 
For an angle 9 between the laser polarization and the 
molecular axis, the effective polarizability can be written 
as a e ff = a± + Aacos 2 #, where a± is the polarizabil- 
ity for the optical field perpendicular to the molecular 
axis and Aa = aii — aj_ is the polarizability anisotropy, 
the difference in polarizability for light polarized parallel 
and perpendicular to the molecular axis. In a linearly 
polarized optical field, the molecules tend to align into 
the field in order to minimize their energy, but because 
of rotational inertia, there is a time lag. With a strong 
optical field, the alignment of molecules causes the linear 
susceptibility to change as a function of time, leading to 
an effective odd-order nonlinearity. To lowest order the 
rotational response is linear in the pulse energy. The case 
of a linearly polarized field has been considered in detail 
previously [21 RB]. 

For a pump and probe of differing polarization, it will 
be shown here that the calculation has additional com- 
plexity. The molecular orientation-dependent dipole mo- 
ment induced by the optical field is p = aE, where a 
is the polarizability tensor. Since the pump and probe 
fields are assumed to be polarized in the xy plane, we 
only need consider the components of the induced dipole 
moment in that plane. These are 



Vx 



Py = 



Aa sin 9 cos 9 cos 0[E P (t)] y 
+ (a ± + Aacos 2 0)[E e (f) + E p (t)] x (11a) 
(a± + Aa sin 2 9) cos 2 (j)[E p (t)] y 
+Aasin6>cos6>cos0[E e (i) +E p (t)] x , (lib) 



where 9 is the angle between the molecular axis and the x 
direction, and </> is the azimuthal angle about the x axis, 
measured with respect to the z direction. For brevity, 
we omit the r argument in this section - note that all 
quantities with a time dependence also depend on r. The 
polarization of the gas is P(i) = N a (p) t , where N a is 
the number density of the gas and ()t denotes the time- 
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dependent ensemble average. We find 



P x = 



Py = 



'- Aa(sin 29 cos i 



[Ep(*)l 



^ N r ,n—ip\,<'J\y 

+N a (a±_ + Aa(cos 2 9) t ) [E e (t) + E p (t)] x 

^ [a x + Aq ((sin 2 6%)] (1 + (cos20) f ) [E p (t)] 
iV 

+ Aa (sin 2(9 cos 4>) t [E e (t ) + E p (t)] x . 



For the calculation we are only interested in the nonlinear 
response, which excludes the static background refractive 
index proportional to N a (a± + Aa/3). Removing this 
static contribution, we find that the nonlinear rotational 
polarization is 
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N a Aa 
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+ 2 

N a Aa 
1 



[E e (*)+E p (t)] x 



sin2^(os<.) t [E p (t)]^ \ , (12a) 
. E p(*)]y 

sin 26» cos 4)t\E e (t) + E p {t)] x 1 . (f 2b) 



We next calculate the time-dependent alignment quan- 
tities (cos 2 6*)t and (sin20cos0)t. The Hamiltonian for 
the interaction of an optical field E and the rotational 
modes of molecules is H = — (l/2)p ■ E. Using Eqs. (11 1, 



we find, to lowest order in the probe field and keeping 
only terms that vary slowly compared with uj e and u> p , 
H '= II ' ' + //" + // ■"'. where 



H xy {t) 



+ [A* e A x e lA]i - r - lAujt + c.c] ) , (13a) 
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- Aa sin 29 cos i 



x [A* e A y e 



iAk-r— iAut 



+ C.C. 



C.C. 
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(13c) 



We seek a perturbative solution for the time evolution of 
the density matrix p using the basis of rotational eigen- 
states \j, m), having quantum numbers j for the total 
rotational angular momentum J and m for the compo- 
nent of J along the x direction. The zeroth order solution 
p(°^ describes the initial thermal equilibrium distribution 
of rotational states. It is diagonal and depends only on 
the total angular momentum j, 



o {0) 

' (jm),(j'm' 



J2 k D k (2k + l)e ^ 



5jj> Sjjijji' , 



(14) 



where Ej — hcBj(j + 1) is the energy of the j rotational 
level (B is a rotational constant, which depends on the 



molecular moment of inertia), and Dj is a degeneracy 
factor related to nuclear spin statistics [14]. For brevity 



We define Pfl = P(jm),(jmy 



The first order solution is 



[H(Oy% m ),Gw)e^' (t '-'W, (15) 



where [] denotes the commutator and Ujj> = (Ej—Ej>)/h. 
We neglect decay of the rotational coherences, as we are 
interested here in the response during the pump pulse 
and a few hundred femtoseconds afterward, whereas the 
rotational coherence lasts > 10 ps at atmospheric pres- 
sure and temperature [Hj. 

Since p^ and arc diagonal, [H, P^'}(j m ),{j'm') = 

\Pjm Pj'm')l- tl \\ \(jm),(j'm') + \Pj m Pj'm' ) U (jrn) ,{j'rn') ' 

where 
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In Appendix B expressions for (jm\ cos 2 9\j'm!) and 
(jm\ sin 26 cos <p\j'm') are given. Finally, we need to cal- 
culate 



(sin26»cos0) t = Tr p {1) {t) sin 26* cos 

P { {jm),{j'm')^M sin 29 cos <j,\ f , m'), 

and the equivalent expression for (cos 2 9) t [H]. In Ap- 
pendix B, it is shown that, in analogy with our previous 
discussion, the ensemble-averaged molecular alignment 
quantities can be split into smooth and grating terms ac- 
cording to (cos 2 0) t = (cos 2 6)1 + 1/3 + ((cos 2 6i)fe iAfe r 



c.c.) and (sin 29 cos 4>)t 
where 
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where 



= 2Aa (j + 1)(j+2) 
3 I5h 



2j + l 
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2j + l 



where pf = (2j + X)pf' As with the plasma response, 



(0) 



the integrals in the grating terms ( 17|18 1 are suppressed 
when Aw is nonzero. However, resonant behavior may 
occur when the pump and probe frequencies differ by the 
spacing between rotational levels. In the experiments 
considered where Aw 7^ 0, Auj is very large compared 
to the spacing between thermally populated rotational 
levels. 

Keeping only terms that lead to a polarization propa- 
gating in the e lkp ' r direction, we find 
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Taking the second derivative with respect to i, we 
find (in the SVEA and keeping only terms without time 
derivatives) 
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-N a Aaujp [(cos 2 e) s t A x (t)e-^ f 



-N a Aauj 2 p [(cos 2 9) s t A y (t)e 

LO 2 
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(19a) 
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The rotational grating produces a probe phase shift that 
approximately follows the pump envelope. In the first- 
order perturbation approximation used here, the rota- 
tional grating contribution is proportional to the pump 
intensity |37) . 



III. CALCULATIONS 

We next perform numerical calculations and compare 
them to previously reported experimental results. Defin- 
ing a local time coordinate r = t — zhq/c, so that 
d/dz — !> d/dz - (n /c)d/dT and d/dt -> d/dr [8], and 
defining barred functions using this local coordinate sys- 
tem A p (r, t) = A p (r, t), Eq. (|2| becomes 



5Ap luJpT _^d 2 P 



NL,p 



Or 2 



(20) 



Note that when these coordinates are used, the probe 
envelope is independent of z in the absence of nonlinear 
interaction. The outgoing probe field is first calculated 
numerically using Eq. (201, with the nonlinear polariza- 
tion sources from Eqs. ( 4|5|10|19 ). The signal measured 
in two types of experiments is then calculated from the 
probe field. 

For the ionization rate w(I) we use the model of 
Popruzhcnko et al. [31], which has been shown to ap- 
proximately agree with TDSE calculations and has been 
used in recent Kramers-Kronig calculations of the non- 
linear response [TH [38] . For the rotational response in N2 
and O2, we use B = 2.0 cm -1 for N2 and B = 1.44 cm" 1 
for O2 and assume T — 293 K. For and Aa we use 
the values reported in |17| . We neglect the vibrational 
response, which in N2 and O2 results in a small (less 
than 10% of the electronic nonlinearity [3S]) effectively 
instantaneous response for the pulse durations used in 
the experiments. 



A. Single-shot spectral interferometry: Degenerate 
vs. nondegenerate probe 

We first simulate the signal in spectral interferome- 
try experiments, which directly measure the time- and 
space-dependent probe phase shift. We calculate using 
the parameters of recent experiments using a chirped 
probe pulse 17, 18J. The pump pulse is assumed to be 
transform limited and Gaussian, centered at 800 nm with 
a full width at half maximum of 40 fs. The probe is ei- 
ther centered at 800 nm (degenerate) or 600 nm (non- 
degenerate). The chirped probe pulse is handled with 
a time-dependent phase in the incident complex probe 
envelope ^4p lc (r, t). For the nondegenerate case, the out- 
going probe field calculated is frequency filtered (as in 
the experiment) to remove components near the pump 
frequency u e before the nonlinear phase shift A$(t) = 
lm{ln{A p (t)/A™ c (t)}} is calculated. 

The time-dependent nonlinear refractive index 
An(t) = n,QUjpA§(t)L/c, where L = 200 /im is the 
interaction length, is shown in Fig. [l] Two different 
calculations are shown: one for N2 at low pump intensity, 
which shows the effect of the rotational response, and 
Ar at high pump intensity, which shows the effect of the 
plasma response. The green dotted line shows the inten- 
sity envelope of the pump pulse, and the thick red (thin 
blue) curves show the response for the probe polarization 
parallel (perpendicular) to the pump polarization. Solid 
and dashed lines show the simulated signal with and 
without the grating response, respectively. The dashed 
line thus indicates the nonlinearity that a laser pulse 
would impart on itself. 

In N2 at intensities below the ionization threshold, the 
bound electronic nonlinearity due to x'' 3 '' leads to a phase 
shift that is proportional to the pump intensity enve- 
lope, and the rotational nonlinearity leads to a delayed 
response peaking about 80 fs after the peak of the pump 
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FIG. 1: Simulated single-shot spectral interferometry signal 
in N2 and Ar for a transform-limited, 40 fs pump pulse cen- 
tered at 800 nm, shown as a green dotted line. The calcu- 
lated nonlinear index An(t) — naCJ p A&(t)L/c, where is 
the time-dependent probe phase shift, is shown assuming a 
nondegenerate probe centered at 600 nm (upper curves) and 
a degenerate probe centered at 800 nm (lower curves). The 
change in index is plotted for the probe polarized parallel 
(thick red) or perpendicular (thin blue) to the pump polar- 
ization. Curves for a nondegenerate probe are offset vertically 
for clarity. The simulated signal is shown including (solid) and 
excluding (dashed) grating effects. 



pulse. When a degenerate probe is used (lower curves), 
an effectively instantaneous rotational grating signal ap- 
pears that, if not properly taken into account in the anal- 
ysis, would make it appear that ni was larger than it is. 
For the probe polarization oriented parallel and perpen- 
dicular to the pump polarization, the rotational grating 
contribution is positive, but not the same magnitude - 
thus, it contributes to a transient birefringence experi- 



ment, as will be shown later. 

In Ar at high intensity, the bound electronic nonlinear- 
ity is accompanied by the nonlincarity due to the gener- 
ation of free electrons (plasma). The "smooth" plasma 
contribution accumulates during the pump pulse and af- 
terward contributes a long-lived, constant phase shift. 
The plasma grating results, for a degenerate probe po- 
larized parallel to the pump polarization, in a strongly 
time dependent negative signal that is a few times larger 
than the "smooth" plasma signal. 

In both cases, when a nondegenerate probe pulse is 
used (upper curves), the only difference between the ex- 
perimentally measured An and the true nonlinear in- 
dex An is the factor of two enhancement in the bound 
electronic contribution to the nonlinearity. These calcu- 
lations provide theoretical justification for the analysis 
performed in recent absolute measurements of the non- 
linearity in the noble gases, N2, O2, and N2O p"7 l 118 ) . 
The extracted values for n2 and Aa are in good agree- 
ment with previous calculations [351 EH EI] 1 experiments 
using less direct methods [39l [42] , and a ratiometric mea- 
surement |43j . Our calculations reinforce the point that 
nondegenerate supercontinuum spectral interferometry is 
the most direct, powerful technique currently available 
for measuring the optical nonlinearity. 

B. Transient birefringence using a degenerate 
probe 

Next we simulate the Loriot et al. degenerate pump- 
probe transient birefringence experiment 144) . A di- 
agram of the heterodyne transient birefringence experi- 
mental apparatus is shown in Fig.[2j The probe is linearly 
polarized at 45° with respect to the pump polarization, 
so A™ c = Ay lc . After the interaction region, the probe 
beam is passed through a phase plate, which provides a 
static retardance 2£ between the x and y probe compo- 
nents, followed by a polarizer oriented at 90° with respect 
to the initial probe beam polarization. The optical field 
after the polarizer is 

A s (Z,t) = e^A x (t)-e^A y (t) = 

e^\A x {t)\e^^ - e -if |i4 1 ,(t)|e < *»<*>. (21) 

In the Loriot et al. experiment, a photomultiplier tube 
is used to detect this field, and the signal is proportional 
to J \A S (£, t)\ 2 dt. The "pure heterodyne" signal is found 
by subtracting the signals for orientations of the phase 
plate producing static retardances £ and — £ [5], 

/OO P OO 

\A s (Z,t)\ 2 dt- / \A S (-H,t)\ 2 dt 
-OO J — OO 

/OO 
\Mt)Uy(t)\ 
-00 

X8in[® x (t)-$ v (t)]dt. (22) 

Since there is little pump-induced change in the probe 
amplitude, \A X \ w \A y \, and Sh is thus proportional 
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FIG. 2: Experimental diagram for transient birefringence ex- 
periment, showing the initial polarization of the pump and 
probe fields, the interaction region inside a gas cell, and the 
detection scheme. 



to the pump-induced birefringence integrated with the 
probe intensity. The heterodyne signal was measured as 
a function of the time delay between the pump and probe 
pulses. We performed 3D numerical calculations of the 
Loriot experiment, assuming that the pump and probe 
beams were Gaussian beams with FWHM waists of 33 
/j,m crossing at 9 = 4°. 



1. Rotational grating effects: simulations for N2 

The calculated transient birefringence signal in N2 and 
O2 at low pump intensity is shown in Fig. |3j with and 
without the rotational grating contribution. The rota- 
tional grating contribution alone is also shown, as well as 
the electronic Kerr contribution alone. Both are roughly 
proportional to the convolved pump intensity envelope. 
The rotational grating contribution in transient birefrin- 
gence is positive, increasing the apparent instantaneous 
nonlinearity observed in a degenerate pump-probe exper- 
iment. 

Table I shows measured values of n 2 from experiments 
by Loriot et al. [5] and from a recent nondegenerate spec- 
tral interferometry experiment by Wahlstrand et al. |17j . 
The reported values for Ar agree within error, but the 
values in N 2 and 2 do not. This can be explained by 
the additional phase shift imparted to the probe by the 
rotational grating in the Loriot et al. experiment (shown 
in Fig. 3]). We used the calculation to subtract the ro- 
tational grating contribution from the Loriot et al. mea- 
surements in N 2 and 2 ; these corrected n 2 values are 
shown in the last column of Table 1. After the correc- 
tion the measurements agree within their uncertainty es- 
timates. 



2. Plasma grating effects: simulations for Ar 

Plots of the calculated heterodyne transient birefrin- 
gence signal versus pump-probe delay in Ar are shown in 
Fig. [4] The signal for 90 fs transform-limited pump and 
probe pulses centered at 800 nm is shown as a function 
of pump intensity in Fig. [4j At low pump intensity the 
induced birefringence is Am — Anj_ — 2n 2 — 2n 2 /3 = 
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FIG. 3: Simulated heterodyne transient birefringence signal 
Sh [Eq. ( 22 1] for N2 and O2 as a function of pump-probe delay, 
assuming transform limited 90 fs pulses centered at 800 nm. 
The peak pump intensity is 10 TW/cm 2 . The curves show 
the calculated signal with and without the rotational grating 
contribution, the contribution from the bound electronic Kerr 
effect alone, and the contribution from the rotational grating 
alone. 



Gas 


Ref. 07] 


n 2 

Ref. rjn] 


(10- 20 cm 2 /W) 

Ref. 10; corrected for rot. grating 


Ar 


9.7 ± 1.2 


10.0 ±0.9 




N 2 


7.4 ±0.9 


11.0 ±2.0 


8.4 ±2.0 


O2 


9.5 ± 1.2 


16.0 ±3.5 


11 ±3.5 



TABLE I: Comparison of ri2 measurements by Wahlstrand et 
al. [XT] and Loriot et al. [9j [TO] . For N2 and O2, values of ni 
from the Loriot et al. experiment are shown with and without 
correcting for the effectively instantaneous contribution from 
the rotational grating. 
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4ri2/3, and the signal is proportional to the pump in- 
tensity envelope convolved with the probe intensity en- 
velope. At high pump intensity, this positive signal due 
to the Kerr birefringence is overwhelmed by the nega- 
tive plasma grating signal. The results are similar to 
those found previously from a one-dimensional calcula- 
tion using a multiphoton ionization rate We find 
that the negative birefringence due to the plasma grat- 
ing appears at a similar intensity to that reported by 
Loriot et al. [U |44j . Note that the normalized intensity 
used by Loriot et al. is the peak intensity divided by 1.7 

0E2]. 

The plasma grating contribution dominates the signal 
at a lower peak intensity than one might expect from re- 
cent measurements using spectral interferometry, where 
the plasma phase shift first appears at approximately 80 
TW/cm 2 [18] . The combination of three effects explain 
this apparent discrepancy. First, the refractive index for 
a given plasma density is larger by a factor of ~ 1.8 for 
an 800 nm probe compared to a 600 nm probe owing 
to the A 2 scaling of the free electron refractive index. 
Thus, while a given plasma density produces a signifi- 
cantly larger phase shift at 800 nm compared to 600 nm, 
the Kerr effect is nearly the same because of its compara- 
tively lesser dispersion. Second, the calculation assumes 
90 fs pulses, about 2.5 times longer than the pulses used 
in |16| I18j . resulting in a higher plasma density for a 
given peak intensity. Finally, the nearly instantaneous 
negative plasma grating signal is a few times larger than 
the "smooth" plasma phase shift at the end of the pulse 

(cf. Fig.g. 

The Loriot et al. data appears roughly symmetric with 
respect to time reversal 04] . Because the plasma grat- 
ing signal is slightly delayed with respect to the Kerr 
signal due to the cumulative nature of the plasma non- 
linearity, the signal is asymmetric at intermediate inten- 
sity. However, note that like two-beam coupling energy 
transfer, the plasma grating phase shift is sensitive to the 
chirp of the pump and probe pulses [25J . If the pulses are 
chirped, the asymmetry of the pump-probe signal is less- 
ened and can even be reversed. This is illustrated in 
Fig. |4j where the calculated signal is shown as a function 
of of group delay dispersion (GDD), assuming that the 
transform-limited pulse width is 90 fs. Both pump and 
probe pulses are assumed to have the same chirp. 



£ 
c 

Z5 

1-2 

TO 
C 

CO 



-6 



1 1 1 1 1 — 

Ar, 90 fs transform-limited 

30 TW/cm 2 — 50 TW/cm 2 



80 TW/cm 2 



2.0 

1.5 

— 1.0 
co 

M — ' 

§ 0.5 

is. 

1 -0.5 

CT> 
CO 

-1.0 
-1.5 
-2.0 




90 TW/cm 2 



100 TW/cm 2 



-200 -100 100 

pump-probe delay (fs) 



200 





Ar, 


80 TW/cm 2 






£\ 700 fs 2 / \ 






\\\ x /rC\\ 




1400 fs 2 


\\ // fs2 






\ y x i800fs 2 






i i i 



-200 -100 100 

pump-probe delay (fs) 



200 



FIG. 4: Sim ulated heterodyne transient birefringence signal 
Sh [Eq. ( 22 1] for Ar as a function of pump-probe delay, assum- 
ing pulses centered at 800 nm. On the left, the calculation 
assumes a Gaussian transform limited 90 fs FWHM pulse, 
and the signal is shown for a few pump peak intensities. On 
the right, the calculation assumes a pump peak intensity of 
80 TW/cm 2 , and the signal is shown for a few values of group 
delay dispersion. Both pump and probe pulses are chirped, 
and the transform-limited pulse width is 90 fs FWHM and 
Gaussian. 



IV. EXPERIMENT 

Based on a numerical calculation in atomic hydrogen, 
Bejot et al. recently argued that the higher-order Kerr 
effect is present only at optical frequencies near the fre- 
quency of an intense driving field |45j , and therefore non- 
degenerate pump-probe experiments, which use a probe 
pulse detuned from the pump pulse frequency, are in- 
capable of observing it. In this section we describe new 
experimental evidence that the negative, nearly instanta- 
neous signal observed using a degenerate probe is caused 



by the plasma grating. 



A. Spectrally-resolved birefringence 

The Loriot et al. experiment measured the temporally- 
convolved transient birefringence as a function of pump- 
probe time delay [9]. Another approach involves spec- 
trally resolving the probe beam. Previously, an experi- 
ment using 400 nm light in air was reported |12) . which 
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FIG. 5: Spectrally-resolved transient birefringence results. The heterodyne-detected signal is plotted as a function of pump- 
probe time-delay and wavelength, (a) Experimental data in Ar for pump pulse energy 42.5 /iJ. (b) Calculation for Ar including 
plasma grating assuming peak intensity 100 TW/cm 2 . (c) Calculation for Ar assuming no plasma grating but a higher-order 
Kerr effect [10] using peak intensity 50 TW/cm 2 . (d) Experimental data in air for pump pulse energy 25 /iJ. (e) Calculation 
for air including plasma grating assuming peak intensity 100 TW/cm 2 . (f) Calculation for air assuming no plasma grating but 
a higher-order Kerr effect [10] using peak intensity 60 TW/cm 2 . 



showed that the power dependence of the negative signal 
was consistent with the plasma grating theory assuming 
multiphoton ionization, which is a good approximation 
at 400 nm, where 4 photons are required for ionization 
of O2. Here, we report new data in Ar and air using 800 
nm pulses. 

The experimental setup is similar to that shown in 
Fig.[2j After the pump and probe interact, a fixed quarter 
wave plate in the probe path produces a static birefrin- 
gence 2£ = 7r/4, and then the probe beam passes through 
a polarizer. The relative polarization of the pump beam 
with respect to the probe is alternated between ±45° to 
make measurements with 2£ effectively 7r/4 and — tt/4. 
The spectrum of the beam after the polarizer is measured 
as a function of pump-probe time delay. Measurements 
using 2£ = tt/4 and 2£ = — 7r/4 are subtracted to gener- 
ate the pure heterodyne spectrally-resolved pump-probe 
signal, with optical frequency on one axis and pump- 
probe time delay on the other. The pump and probe 
pulses were 42 fs FWHM, approximately transform lim- 
ited pulses. Results are shown in Fig. [5^, for Ar and 
Fig[5tl for air. 



To simulate the experimental data, we numerically cal- 
culate the optical field at the output of the polarizer, 



given by Eq. (21 1, and then Fourier transform it and 
calculate the magnitude squared. The pure heterodyne 
spectrum is found by subtracting the calculated Fourier 
transform for 2£ = ir/4 and 2£ = — 7r/4. The signal 
is plotted as a function of optical frequency and pump- 
probe time delay in Fig. bjp (Ar) and Fig. [5^ (air) 
for the plasma grating model and in Fig. [5p (Ar) and 
Fig. [5^ (air) neglecting ionization but including higher- 
order Kerr coefficients [ID]. The shape of the signal 
clearly matches the plasma grating model better than the 
higher-order Kerr effect model. The difference between 
Fig. [5]d and [5]: can be attributed mostly to the cumula- 
tive nature of the plasma nonlinearity, which also results 
in the slight asymmetry in the spectrally integrated signal 
at high intensity observed in Fig. [2^,. This asymmetry is 
more pronounced in spectrally resolved data than in the 
pump-probe traces from the Loriot et al. technique [5J. 
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FIG. 6: Degenerate spectral interferometry results, (a) Ex- 
perimental setup, (b) Extracted time domain phase shift as 
a function of the pump pulse energy for parallel polarizations 
(solid) and perpendicular (dashed). The curves for different 
pump pulse energies are vertically offset for clarity. 



B. Degenerate, chirped pulse spectral 
interferometry 

To further investigate the origin of the negative bire- 
fringence, we use spectral interferometry with a nearly 
degenerate probe pulse. The principle of the degenerate 
experiment is identical to our previous nondegenerate im- 
plementation using a supercontinuum probe [16] except 
that the pump and probe are not collinear. This allows 
us to orient the pump polarization parallel to the probe. 
A sketch of the experimental apparatus is shown in Fig. 
[6ji. The pump is a 40 fs FWHM pulse centered at 800 
nm. The pump and probe beams cross at 3°, allowing the 
pump beam to be blocked before probe detection. The 
two beams are crossed inside a vacuum chamber back- 
filled with 0.5 bar of Ar. The probe and reference pulses 
are chirped using a block of glass to a group delay dis- 
persion of 1650 fs 2 . 



Results of the experiment are shown in Fig. [6Jd. For 
the pump beam polarized perpendicular to the probe 
beam, the result is identical to that found with collinear 
geometry [IB] . We observe a positive instantaneous phase 
shift near the time of peak pump intensity (t — 0) due to 
the optical Kerr effect, and a negative, long-lasting signal 
from electrons freed by ionization. When the pump beam 
is polarized parallel to the probe beam, we in addition 
observe a negative signal near t = that dominates the 
positive Kerr signal at high intensity. We attribute the 
oscillatory signal to interference between supercontinuum 
generated by the pump pulse and the probe pulse. It is 
important to note that as the intensity is increased, the 
negative signal at t — appears simultaneously with the 
smooth plasma signal. This is strong evidence for the 
plasma grating interpretation. In addition, the shape of 
the signal is consistent with our calculation, shown in 



V. CONCLUSION 



We have presented a detailed theory and numerical cal- 
culations of pump-induced phase shifts of probe pulses in 
media with electronic, plasma, and rotational nonlineari- 
ties. Such experiments are typically used to measure the 
nonlinear response of a material. The theory includes the 
effect of the nonlinear interference grating formed by the 
space and time overlap of the pump and probe fields. We 
showed that proper interpretation of the probe response, 
in terms of a medium's fundamental nonlinearities, can 
only be achieved by correctly accounting for the effect of 
the interference grating. In particular, our calculations 
reveal how the presence of nonlinear interference gratings 
(in both the plasma and rotational responses) in recent 
pump-probe experiments have resulted in misinterpreta- 
tion of the obtained results for the nonlinear response of 
gases at high intensity [S]. 

In addition, we have presented results from two new 
degenerate pump-probe experiments, spectrally resolved 
transient birefringence and noncollinear spectral inter- 
ferometry, for the specific purpose of investigating the 
effect on a probe pulse of the nonlinear interference grat- 
ing. The results of both experiments are well reproduced 
by our calculations. The results of this paper reinforce 
the idea that, for the purpose of measuring the nonlinear 
response of a medium using cross phase modulation, a 
nondegenerate pump-probe experiment is preferable to a 
degenerate one. As shown earlier, for a nondegenerate 
experiment the only difference in the measured nonlinear 
refractive index and the nonlinear refractive index in- 
duced by a pulse and itself is a factor of 2 enhancement 
in the instantaneous, bound electronic nonlinearity. 
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Appendix A: Electronic Kerr effect 

Putting the total optical field E = E p + E e into Eq. ^ and recalling that the pump is polarized along x and the 
probe is polarized in the xy plane, we find (keeping only terms containing g^p-""-'"? 4 ) 



(P e i)x(r,t) 

(Pel)yM) 



27T 

n c 
2tt 
n c 



^XxxxxAxi? , t) + (yXxyxx Xxxyx Xxxxyj AyijC, ^e( r i t) 
3XyxxxAx( r , t) + \Xyyxx Xyxyx + Xyxxyj A y {v , t) I e (r, t)e 



ikn-r— iuj v t 



(3) (3) 

where we have rewritten the expressions in terms of the pump intensity. In an isotropic medium, \xyxx — Xxxyx = 

(3) (3) 

Xxxxy = Xyxxx = 0, so we neglect all terms containing those elements. Taking the second derivative with respect to 
time, we find 



4TTd 2 (P cl ) x _ 2n n 2 d 2 



dt 2 



c 2 Q t 1 



(leA 



where n<i is defined in the text. Expanding the derivative, we find 



dt 2 



dt dt 



dt 



For many-cycle pulses, only a few_terms are important above. In the SVEA, second-order and higher time derivatives 
are neglected. This lead 
similarly derive Eq. ([5]). 



are neglected. This leads to Eq. (|4j). In an isotropic medium, Xmxx — Xyxyx — X^yxxy — X^xxxx/i, from which one can 



Appendix B: Rotational grating 



The quantity (j,m\ sin 28 cos <p\f , to') can be expressed in terms of spherical harmonics, 

(J, m\sm26 cos <p\ f , to') = (j,m\ sm#cos6>(e^ + e~^)\j' ,m') 

F8w 



({j,m\Yj(e,4>)\j',m') - VMYiHWUW)) 



These integrals of 3 spherical harmonics are [33] 



(j,m|sin20cos0|j , ) m') = {[A$(j')8, 



m,m' — l -f- 



A + m (j')5 m , m > + i I 5 



3 J' +2 



+ [AFV')8 m , m ,-i - Aa m (j')8 m , m , +1 ] S jd > + [A m (f)S m ^^ - AZ m (j')8 m , m >+i] hi>-*} > ( B1 ) 



where 



A'm 

Ah n (3) 
A m (j) 



(j - m)(j - to + l)[(j + 2) 2 - m 2n 1/2 



(2j + l)(2j + 3) 2 (2j+5) 
(1 + 2m) [(j -m)(j + l + m)] 1 / 2 
(2j - l)(2j + 3) 

(j + m + l)(j + m)[(j - l) 2 - to 211 1/2 
(2j-3)(2j-l)2(2j + l) 



Similarly, 



(j,m| cos 2 9\f,m') = B™(f)5 jtr+2 5 mm , + B™(f)5 jd ,- 2 5 mm , + 



W) 



(B2a) 
(B2b) 
(B2c) 

(B3) 
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where 



B?(j) = 
B n l (j) = 
B m (j) = 



[(j+2) 2 -to 2 ][(.? + 1) 2 -to 2 



-,1/2 



(2j + l)(2. ? + 3) 2 (2j + 5) 

2 j(j + 1) - 3to 2 
3(2j-1)(2j+3)' 

Q- 2 -m 2 )[(j-l) 2 -m 2 ] i 1/2 
(2i-l)(2i + l)2(2j-3). 



(B4a) 
(B4b) 

(B4c) 



Using Eqs. ( |Bl|B3[ ) in Eqs. ( |l3b|l3c| and those plus Eq. §U§ in Eq. ([TH]), we find that only certain terms in 



are nonzero. These are 
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p 0',m),(j+2,m-i) 



(i-2)( 

-B+Ci + 2) (p) 
-^ r + (J - 2) (p 

0'"2) (, 
-A m (j + 2) (p< 
A = m 0" + 2) (, 



.4 



,(o) 
i 

(0) 
3 



(0) 



(o) 



(0) 

Pi 



pf 2 )G^- 2 (M) 
pf 2 )G^_ 2 (r,t) 
Pj+2) Gj,j+2(r,t) 
pf 2) G jt3+2 (r,t), 



where 



~2?T 

iAa 
~4h 



e i UjJ , (t'-t) [\ Ae \ 2 (r,t') + (A* e (r,t')A x (r,t')e' 



Ak-r-iAwt' 



+ c.c.)]d£. 



Finally, 



(sin26»cos0) t = Tr(p sin 29 cos <j>) = 53 p[V m) (f m , } (r, t)(j, m| sin26> cos <j)\ f , to') 

= E [-^= m "Hi)p[il ) , ( ,_2, ro+1 )(r,i)+^- 1 (i)p!il ))(i _ 2jm _ 1) (r,t) 



-A] 



0')^(j,m),(j'+2,TO+l)( r ' + ^+ i (j')^(j,m),(i+2,m-l)( r '*) 

C?')P(j,m),0'-2,m+l)( r '*) + ^- 1 (j')' J (j,m),(i-2,m-l)( r '*) 



(1) 



2 [-a:™- 1 

_ ^+ m (i)/ ? (j,m-l),(i+2,m) + ^■+0')/ ? ( i3 ',m+l),(j+2,m)( r '*) 



Inserting Eqs. (B5I into this, we find 

(sin20cos0) t = ^ [AZ m -\j)A™{j - 2) (pf - pf 2 ) G JJ _ 2 (r,t) 

I , m 

- 2) (pf > - pfl 2 ) G w _ a (r,t) 
- 2)A m (j + 2) (pf - pf 2 ) CWr.t) 
- 2)Al" l (j + 2) (pf pf 2 ) G j+2J (r,t) 
= E [^ m_1 C?)^+(j - 2) + ^(fl^Ttf - 2)] {pf pf 2 )i^_ 2 (r, t) 



where in the last line we have used 



G jr {r,t)-Gj,j{T,t) =H jf (r,t) = 



Aa 

~2h 



sm[uj jd ,(l/ - t)][A* e (r,t)A y (r,t)e 



i Akr- iAujt' 



■]dt'. 



(B5a) 
(B5b) 
(B5c) 
(B5d) 
(B5e) 
(B5f) 



(B6) 



14 



Using Eqs. (B2| in Eq. (B6| and cancelling many factors 



leads to Eq. (18). A similar calculation, which we do not 



reproduce here because it has been described previously 
QI], leads to Eqs. ( pjl7| . 
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